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Ground State and Magnetization Process of the Mixture of Bond- Alternating and 
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The mixture of bond-alternating and uniform S = 1/2 antiferromagnetic Heisenberg chains 
is investigated by the density matrix renormalization group method. The ground state mag- 
netization curve is calculated and the exchange parameters are determined by fitting to the 
experimentally measured magnetization curve of C\iC\2 X Bt2(i-x) (7-pic)2 • The low field be- 
havior of the magnetization curve and low temperature behavior of the magnetic susceptibility 
are found to be sensitive to whether the bond-alternation pattern (parity) is fixed all over the 
sample or randomly distributed. The both quantities are compatible with the numerical results 
for the random parity model. 
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1. Introduction 

In recent studies of quantum many body problem, the 
ground state properties of the random quantum systems 
have been attracting a renewed interest. Among them, 
the effect of randomness in one-dimensional quantum 
spin systems has been extensively studied theoretically 
and experimentally. 1-13 ' 

The real space renormalization group (RSRG) analysis 
of the S = 1/2 random bond antiferromagnetic Heisen- 
berg chain 1-3 ) has shown that the ground state of the 
5=1/2 uniform antiferromagnetic Heisenberg chain 
is unstable against the bond randomness of infinitesi- 
mal strength and the random singlet phase is realized 
in which the spins form singlet pairs randomly with dis- 
tant partners. 3 ' In this phase, the spin-spin correlation 
decays by the power law with exponent —2 and log- 
averaged energy gap scales as < InA >~ —a^/N +const., 
where N is the system size and a is a constant. This has 
been confirmed numerically 4 ' by means of the density 
matrix renormalization group (DMRG) method. 14 ' 

On the other hand, the effect of randomness on the 
spin gap state of one-dimensional quantum spin systems 
has been extensively studied 5-11 ' related with the dop- 
ing effect in bond-alternating chains 8 ' 9 ' and spin lad- 
ders. 10 ' Hyman and coworkers 5,6 ' have applied the 
RSRG method to the spin-1/2 random dimerized an- 
tiferromagnetic Heisenberg chain and have shown that 
the dimerization is a releveant perturbation to the ran- 
dom singlet phase. 5 ' 6 ' This implies that the ground 
state of the spin-1/2 dimerized antiferromagnetic Heisen- 
berg chain is stable against infinitesimal randomness. 
They also argued that the ground state of this model 
is the random dimer phase which belongs to the quan- 
tum Griffiths phase in which the spin-spin correlation 
length is finite and log-averaged energy gap scales as 
< InA >^ — zlniV + const, with finite dynamical ex- 
ponent z. This is also verified by the DMRG calculation 
by the present author. 7 ' 

Motivated by the quite distinct nature of the uniform 



chain and bond-alternating chain against randomness, 
Ajiro and coworkers experimentally investigated the 
magnetic properties of the compound CuCl2 a; Br2(i_ 2 ;)(7- 
pic)2 12 ' 13 ' which is the mixed compound of S = 1/2 
uniform and bond-alternating antiferromagnetic Heisen- 
berg chains. The x = 1 compound CuCi2(7-pic)2 is an 
S =1/2 bond- alternating antiferromagnetic Heisenberg 
chain 15 ' and the x — compound CuBr2(7-pic)2 is an 
5 = 1/2 uniform antiferromagnetic Heisenberg chain. 12 ' 
In order to understand the low temperature magnetic 
properties of these compounds, we theoretically investi- 
gate the ground state and low energy properties of the 
mixture of bond-alternating and uniform 5 = 1/2 an- 
tiferromagnetic Heisenberg chains using DMRG method 
in the present work. 

This paper is organized as follows. In the next section, 
the thoretical model which describes the mixed com- 
pound CuCl22;Br2(i_ x )(7-pic)2 is presented. In section 3, 
the magnetization curve is calculated and the exchange 
parameters are determined from overall behavior of the 
magnetization curve. The numerical results of the low 
energy magnetic excitation spectrum is presented in sec- 
tion 4. Based on the finite size scaling analysis of the 
log-averaged energy gap, the low temperature singular- 
ity of the magnetic suscsptibility is predicted and the 
explanation of the experimentally observed temperature 
dependence of susceptibility data is given. The last sec- 
tion is devoted to summary and discussion. 

2. Theoretical Model of CuC^B^i— x ) (7-pic)2 

In the compound CuCl2xBr 2 (i_ 2; )(7-pic)2 , the Cu-Cu 
bond is bibridged by CI and/or Br ions, so that there ex- 
ist Cu <^}> Cu , Cu <q> Cu and Cu <g£> Cu bonds 

in this compound. In the pure compound CuCl2(7-pic)2 
PI 

, two kinds of Cu <q> Cu bonds alternate along the 
chain, although this material appears to be a uniform 
chain from chemical formula. This bond alternation is 
induced by the freezing transition of rotational motion 
of methyl-group at 50K. 15 ' For finite Br-concentration 
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(0 < x < 1), this compound can be regarded as an as- 
sembly of the finite length clusters connected only by the 
Cu <q> Cu bonds (hereafter called 'bond-alternating 

cluster'; abbrivaiated as BAC) with Cu <g£> Cu and/or 

Cu <q> Cu bonds in between. However, it is not ob- 
vious whether the rotational order of methyl-group re- 
mains long ranged even in the mixed compound. If this 
order is long ranged, the bond alternation pattern is com- 
mon among different BAC's even though they are sep- 
arated by Cu <g£> Cu or Cu <q> Cu bonds. In 
what follows, this case is called fixed parity case. On the 
other hand, if the rotational order of methyl-group is 
also cut into short range order by Br-substitution on CI 
sites, the bond alternation patterns of the Cu <q> Cu 
bonds are uncorrelated among different BAC's once they 
are separated by Cu <q> Cu or Cu <gj.> Cu bonds. 
This case is called random parity case. For the moment, 
there is no experimental evidence which case is realized 
in CuCl2a;Br2(i_ 2 ;)(7-pic)2 . In the following, we there- 
fore investigate both cases theoretically. 

The strength of these bonds are denoted by J (strong 
Cu <q> Cu bond), aJ (weak Cu <q> Cu bond), 

J' (Cu <Br> Cu bond) and J" (Cu <|j> Cu bond), 
respectively in the following Hamiltonian which describes 
the mixed chain. 



Regular bond alternating chain 



in phase 



in phase 



Fixed Random 
Parity Parity 
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0% 50% 
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Fig. 1. Schematic view of the local bond configuration. The thick 
lines denote the Cu <qj> Cu or Cu <g£> Cu bonds and the 
thin single and double lines, the alternating Cu <(j}> Cu bonds. 
The uppermost configuration R corrsponds to the regular bond- 
alternating chain (CuCl2(7-pic)2 )• Two possible local bond 
configurations in the mixed chain CuCl2^Br 2 (i_ a; ) (7-pic)2 are 
shown as A and B. The clusters in the rectangles arc bond- 
alternating clusters (BAC's). In the random parity case, the 
type- A and type-B configurations occur with the same weight. 
In the fixed parity case, only type-A configuration occurs. 



3. Magnetizaton Process 

From the measurement for the pure systems (x = 
0, 1), the coupling constants J, a and J" are determined 

as, 12 - 13 ) 



H 



N 

£ 



2JiSi Si+i, 



(1) 



J = 13. 2K, J" = 20.3K, a = 0.6. 



(3) 



Ji = < 



J 

aJ 
J' 
J" 



for i + ip 
for i + ip 



even 
odd 



P = x , 

p = 2x(l — x), 
p=(l-x) 2 



(2) 

where p is the probability of realization of each type of 
bonds. We define the integer 'parity' i p for each BAC. 
Corresponding to the fixed and random parity cases ex- 
plained above, the parity i p takes the following values in 
each cases. 

(1) Random parity case : i p takes the values or 1 
randomly in different BAC's. 

(2) Fixed parity case : i p = for all BAC's. 

The possible local bond configurations realized in each 
case can be visually explained using Fig. 1. The upper- 
most configuration R shows the regular bond-alternating 
chain corresponding to CuCl2(7-pic)2 . The local config- 
urations of type-A and B can occur in the mixed chain 
CuCl2 2: Br2(i_ 2: )(7-pic)2 . The two BAC's in type-A con- 
figuration are in phase with regular chain R. On the other 
hand, in type-B configuration, the left BAC in the solid 
rectangle is in phase with R while the right BAC in the 
dotted rectangle is out of phase from R. The local con- 
figurations of type-A and B (with its mirror inversion) 
occur with the same weight in the random parity case 
but only type-A configuration occurs in the fixed parity 
case. 



Using these parameters, we have calculated the magne- 
tization curve of the mixed chain using DMRG method 
for various choices of J' and determined J' / J = 1.3 
so as to reproduce the overall magnetization curve of 
CuCl2a;Br2(i_ 2 ;)(7-pic)2 . The chain length N is fixed 
to 100 and the average is taken over 100 samples. The 
maximum number of the states kept in each subsystem 
in the course of DMRG calculation is 300. It should be 
noted that the overall magnetization curve does not de- 
pend whether the parity is fixed or random, except for 
the low field part where the magnetization curve is more 
convex in the random parity model as shown in Fig. 2. 
The low field part of the magnetization curve reported in 
refs. 12 and 13 seems to be consistent with the random 
parity model. However, due to the finite temperature ef- 
fect in the experimental data and finite size effect in the 
numerical calculation, the direct quantitative compari- 
son is difficult. Actually, the most significant difference 
between the two models is the low field singularity of the 
magnetization determined by the low energy singularity 
of the magnetic excitation spectrum. At finite tempera- 
tures, however, this singularity is rounded to yield finite 
magnetic susceptibility which diverges or vanishes as the 
temperature is lowered. Therefore, for the quantitative 
comparison with experiment, we concentrate on the low 
temperature singularity of the susceptibility which can 
be deduced from the finite size scaling analysis of the low 
energy magnetic excitation spectrum in the next section. 
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Fig. 2. Magnetization curves of the fixed parity model (solid line) 
and random parity model (dotted line) with x = 0.6. The step- 
wise lines are the magnetization processes of pure systems (x = 
and x = 1). For x = 0.6, the error bars are comparable with the 
width of the lines. 



4. Low energy magnetic excitation spectrum 

4--1 Fixed parity case 

The lowest excitation gaps with total magnetization 
gtot _ 2 f or eacn sample is calculated using the infinite 
size DMRG method 4 ' for 10 < N < 80 (x < 0.5) and for 
10 < N < 160 (x > 0.6). The average is taken over 200 
samples. The number of the states kept in each subsys- 
tem during the course of the DMRG calculation is 60. 
The accuracy is confirmed by checking the coincidence 
of the energy gap with (<Stot, £>tot) = (1) 0) and that with 
(Stot, <Stot) = (Ijl) taking into account the SU(2) sym- 
metry of the present system. 

The log-averaged energy gap < InA > is plotted 
against hiN in Fig. 3 which show a fairly good linear 
behavior 

< InA >= lnAo - zlnN, (4) 

for large N and x > 0.4. For x < 0.3, the data are 
not well fitted by (4). Presumably, in this regime, the 
concentration of the alternating bond x 2 is too small 
and true asymptotic behavior is not observable for the 
present system sizes. This result implies that the ground 
state of this model is the quantum Griffiths phase char- 
acterized by the finite dynamical exponent z at least for 
x > 0.4. In the quantum Griffiths phase, the low temper- 
ature susceptibility x{T) and low field magnetization at 
zero temperature M(H) should behave as x(T) ~ j^h—i 
and M{H) ~ H 1 / 2 . 5 '^ 

Physically, this result is reasonable. In this model, the 
dimerization pattern is fixed over the whole system. As 
pointed out by Hyman and coworkers, 5 ' the dimerization 
is a relevant perturbation to the random singlet phase 
and drives the system into the random dimer phase which 
belongs to the quantum Griffiths phase. 



In the experiment, it is found that the low temper- 
ature susceptibility of CuCi2;cB r 2(i— a;) (7 _ pi c )2 shows a 
divergent behavior as \ ~ y/ 3-1 with an exponent 
(3 ~ 0.67. 12 - 13 ' 16 ' The value of (3 slightly depends on 
the concentration x. At first sight, the above result ap- 
pears to be consistent with the prediction of the fixed 
parity model, if we identify z = 1/(3 as proposed in ref. 
12. Based on this idea, we determine the exponent z by 
fitting < InA > in Fig. 3 by Eq. (4) . The result is shown 
in Fig. 4 in terms of (3 = 1/z. Although the value of (3 
for small x is difficult to determine from the numerical 
data, we may expect it would roughly behave as shown 
by the dotted line considering that (3 = 1 for x = 
(Luttinger liquid). The exponent f3 is, however, always 
larger than or around unity as shown in Fig. 4 within 
the regime x > 0.4. Therefore the fixed parity model 
does not explain the experimentally observed exponent 
(3 ~ 0.67. Actually, if we look into the experimental 
data closely, the exponent (3 turns out to decrease with 
x, although the temperature of the experiment is not low 
enough to exclude the contribution of higher energy ex- 
citations. On the contrary, in the fixed parity model, (3 
increases rapidly with x for x > 0.4 as shown in Fig. 
4. The physical reason of the sharp increase of /3 is ob- 
vious. In the limit x — > 1, the system is gapped, so that 
z should tend to zero. 
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Fig. 3. System size dpendence of < InA > in the fixed parity case. 
The gap is measured in units of J. 



4-2 Random parity case 

In this case also, the lowest excitation gaps with total 
magnetization 5* ot = 1 for each sample is calculated 
using the DMRG method for TV < 80 in the same way 
as in the preceding subsection. The maximum number 
m of the state kept in the course of DMRG calculation 
is m = 200. 

The log-averaged magnetic energy gap < InA > is 
plotted against \A/V in Fig. 5 which show a fairly good 
linear behavior 

< InA >= lnA - aVN. (5) 
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"0 0.5 x 1 

Fig. 4. The exponent fi plotted against x in the fixed parity case. 



This implies that the ground state is the random sin- 
glet state in which the low temperature susceptibil- 
ity x{T) and low field magnetization at zero tempera- 
ture M(H) should behave as x(T) - l/(T(lnT) 2 ) and 
M(H) ~ (\n(l/H))- 2 .V In this case, the dimerization 
pattern is not long ranged. Therefore the translational 
symmetry is preserved on average and the random singlet 
phase remains stable. 

However, it should be noted that this asymptotic form 
for x is oni y valid in the true low temperature limit. 
At finite temperatures, we may define the temperature 
dependent effective value of z even in the random parity 
case in the following way. 
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Fig. 5. System size dependence of < InA > in the random parity 
case. The gap is measured in units of J. 

Let us define the size dependent effective value of z by 



ZeS 



d < InA > 
rflniV ' 



(6) 



based on the formula (4). This size dependent exponent 
can be interpreted as the energy dependent effective ex- 
ponent because the typical energy scale of the cluster is 
related with the cluster size N by the relation (5). Thus 
we find that z e g is given by 



N 



Therefore, in the temperature range around T, we can 
define the temperature dependent effective value of z by 



(8) 



by setting < InA >= ln(fc B T) in (7). 

From the fit to (5) we have determined the value of 
Ao. Around the temperature range where the measure- 
ment is made (T ^> 2K), the values of the effective 
exponent (3 c g(T) = l/z e g(T) are plotted against T in 
Fig. 6 for avrious values of x. They are also plotted 
against x for T = 2K and 4K in Fig. 7. It is clearly 
seen that (3 e g ranges from 0.4 to 0.7 for 0.8 > x > 0.2 . 
This explains why the measured susceptibility exponent 
is around these values. Especially, at the fixed tempera- 
ture, the exponent (3 decreases with x in accordance with 
the experimental tendency. 

Physically, the rr-dependence of f3 e g is understood as 
follows. The quantity Ao is essentially the bare energy 
scale of the system. In the spirit of RSRG method for 
random quantum spin chains, 1 ~ 3 - ) it is of the order of the 
largest exchange coupling 2 J" <~ 40. 6K which is deci- 
mated in the first step renormalization. Therefore we 
may roughly estimate as (3 e g ~ 2/ln(2 J" /UbT) ~ 0.66 
at T = 2K and this value is insensitive to x. However, as 
x appraoches zero, the ground state should tend to the 
Luttingcr liquid characterized by z = f3 = 1. Therefore 
/3 e ff increases to approach this value with the decrease 
of x. Actually, for small x, there exist many large uni- 
form clusters coupled with the strongest bond J" so that 
the first step renormalization cannot start with the sin- 
gle bond. Therefore, the effective value of A is reduced 
from J" and @ e g is enhanced. 




Fig. 6. The effective exponent /3 e ff plotted against temperature 
TK. 



5. Summary 

The magnetization process of the mixture of uniform 
and alternating-bond S = 1/2 antiferromagnetic Heisen- 
berg chain is calculated using the DMRG method in the 
ground state. The magnetization processes of random 
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Fig. 7. The effective exponent f3 e g plotted against CI concentra- 
tion x. 

and fixed parity cases are similar to each other except 
for the low field singularity which reflects the low energy 
spectrum. Comparing with the experimental magnetiza- 
tion curve for CuCl 2x Br2(i_ x )(7-pic)2 the exchange pa- 
rameters are determined. The low field magnetization 
curve of CuCl2 X Br 2 (i_ x )(7-pic)2 looks similar to that of 
the random parity case although the direct comparison 
is difficult due to the finite temperature effect in experi- 
ment and finite size effect in numerical calculation. 

For the quantitative analysis of the low field limit, 
the low energy energy spectrum obtained by the DMRG 
method is analyzed using the finite size scaling method. 
The ground state is the random singlet phase for the 
random parity model and it is the quantum Griffiths 
phase in the fixed parity model. The low temperature 
singularity of the experimentally observed susceptibility 
of CuCl2 2: Br2(i_ a ;)(7-pic)2 is consistent with the theo- 
retical results for the random parity model assuming the 
temperature dependent effective exponent f3 c g. 

It should be noted that our calculation predicts that 
the effective value of (3 decreases with the CI concentra- 
tion for a; < 0.5. Although this seems to be compatible 
with the experimental data, more precise measurement 
of susceptibility at lower temperatures is hoped in the 



future including the possibility of direct measurement of 
the random singlet behavior x(T) ~ l/(T(lnT) 2 ). 
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